Graphene on metal surface: gap opening and n-doping 
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Graphene grown on metal surface, Cu(lll), with a boron nitride(BN) buffer layer is studied for 
the first time. Our first-principles calculations reveal that charge is transferred from the copper 
substrate to graphene through the BN buffer layer which results in a n-doped graphene in the 
absence of a gate voltage. More importantly, a gap of 0.2 eV which is comparable to that of a 
typical narrow gap semicondutor opens just 0.5 eV below the Fermi-level at the Dirac point. The 
Fermi-level can be easily shifted inside this gap to make graphene a semiconductor which is crucial 
for graphene-based electronic devices. A graphene based p-n junction can be realized with graphene 
eptaxially grown on metal surface. 
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Graphene, first isolated by Novoselo^, have attracted 
a lot of research interest recently because of its intriguing 
physics for fundamental studies and its potential appli- 
cations for the next generation electronic devices, such as 
novel sensors 2 - and post-silicon electronics^^ As a 2D 
crystal, graphene was found not only to be continuous 
but to exhibit high crystal qualityjii^ in which charge 
carriers can travel thousands of interatomic distances 
without being scattered] 9 ' 10 ! 11 In addition, due to its spe- 
cial honeycomb structure, band structure of graphene ex- 
hibits two intersecting bands at two inequivalent K points 
in the reciprocal space and its low energy excitations are 
mass-less Dirac fermions near these K points because of 
its linear (photon-like) energy-momentum dispersion re- 
lationship. This results in very high electron mobility 
in graphene which can be further improved significantly, 
even up to w 10 5 cm 2 /V-s. It also allows controls of 
carrier type (electron-like or hole-like) and density by 
electric-field^. These are different from conventional dop- 
ing of semiconductors, for example via ion implantation. 
With such features, graphene is promising in many ap- 
plications, such as bipolar devices which comprises junc- 
tions between hole-like and electron-like regions, and p-n 
junctions which can be configured by gate- voltage within 
a single atomic layer . 12 ' 13 ' 14 Graphene based devices can 
be expected to have much more advantages than silicon- 
based devices. 

At present, there are mainly two methods of produc- 
ing graphene samples. In the first, an almost freestand- 
ing graphene is produced by mechanically splitting off 
bulk graphite crystals and depositing it onto a SiO^/Si 
substrate^ Although this method is very convenient, it 
is difficult to produce high quality graphene structures. 
In the second method, an ultrathin graphene layer is 
formed by vacuum graphitization through silicon dcplc- 
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tion from a SiC surfaced The large lattice mismatch 
between SiC substrate and graphene results in graphene 
produced using this method being fragile and contain- 
ing many defects.— It is thus desirable to find a way 
of fabricating graphene which can overcome the above 
disadvantages. Pristine graphene is a zero-gap semicon- 
ductor because its Fermi level exactly crosses the Dirac 
point. An energy gap is required for many practical ap- 
plications. 

A gap can be induced in graphene by a weak interac- 
tion with the substrate when it is grown on a suitable 
substrate.— However, choosing the right substrate is not 
trivial. If the interaction is strong, the advantages of 
graphene will be lost which is the case when a single 
layer of graphene grown on SiC— On the other hand, 
a too weak interaction between graphene and the sub- 
strate is not able to change the property of graphene and 
to open a gap^ Attempts were made to grow graphene 
directly on metal substrate , 21 ! 22 which provides an inter- 
action stronger than van der Waals interaction (such as 
graphene grown on h-BN substrate^) but weaker than 
covalent interaction (such as that between graphene and 
a SiC substrate^ 9 -), which indeed led to opening of band 
gap. In addition, large area of flat single-crystal metal 
surface is easy to obtain experimentally. Therefore, metal 
substrate can be useful in graphene-based electronic- 
devices. The only drawback is that the energy bands 
and the Dirac cone of graphene are submerged in the elec- 
tronic sea of metal and the unique property of graphene 
cannot be utilized when it is grown directly on metal 
substrate. It is thus important to find proper substrate 
to overcome this disadvantage. In this work, we investi- 
gate graphene grown on Cu(lll) substrate with a boron 
nitride (BN) buffer layer between them. The experimen- 
tal lattice constant of Cu(lll) surface is 2.55 A which is 
very close to those of the BN sheet (2.50 A) and graphene 
(2.46A). The slight mismatch is acceptable since the C-C 
bond and B-N bond are very strong. It is noted that in 
the growth of graphene directly on SiC, the excess car- 
bon atoms self-organize to form the honeycomb structure 
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which overwhelms the covalent bonding with substrate, 
even through the lattice mismatch between the graphene 
layer and the SiC substrate is as high as 8%.— Thin car- 
bon layer consisting of sp 2 -hybridized carbon modifica- 
tion in the form of graphene stacks over copper nanopar- 
ticles have been produced.— Heteroepitaxial graphene/h- 
BN double layer has also been realized on some metal 
substrate by chemical vapor deposition(CVD)^ It can 
be expected that including a BN buffer layer will im- 
prove the stability of the heteroepitaxial graphene/h-BN 
structure on Cu(lll) surface. Furthermore, the special 
electronic structure of this system ensures that it will be 
more suitable for applications in electronic devices, which 
will be further discussed below. 

First-principles calculations based on the den- 
sity functional theory were performed using the 
VASP code^^i^L^ The projector augmented wave 
pscudopotentials22iwere used for electron-ion interactions 
while the local density approximation (LDA) was used 
for exchange-correlation function. A special 36 x 36 x 1 k- 
point sampling was used for the surface Brillouin-zone 
integration. The plane wave basis set was restricted by 
a cutoff energy of 400 eV. The in-plane lattice constant 
of the heteroepitaxial graphcnc/h-BN/Cu(lll) structure 
was set to 2.49 A which is the lattice parameter of 
Cu(lll) determined from our first-principles calculation 
within LDA. In structural relaxation, five atomic layers 
of Cu were used to model the substrate, with atoms in 
the bottom two layers fixed to their bulk positions and h- 
BN and graphene adsorbed to the top surface. Electronic 
structure calculations were carried out using a ten layer 
substrate model with two center layers fixed. Graphene 
and BN layer are adsorbed on both sides. A vacuum 
region of at least 18 A was used to avoid interaction be- 
tween the top and bottom surfaces. All structures were 
fully relaxed, and the forces acting on an atom is < 0.03 
eV/A in the relaxed structures. 

The graphene and BN layers were placed on the lxl 
unit cell of Cu(lll). There can be different atomic ar- 
rangements for BN and graphene relative to the sub- 
strate. First of all, the BN sheet can have six different 
arrangements over Cu(lll), with either B or N on the 
fee or hep site, and atop a Cu atom of the substrate, 
respectively. Among these, two arrangements, with N 
on the atop site while B on either the fee or hep site as 
shown in Fig. 1(a) and Fig. 1(b), respectively (referred as 
structures Al and A2 respectively in the following), are 
energetically favored. They are very similar to the struc- 
ture of BN sheet on Ni(lll).^ The energies of these two 
structures are almost the same and are lower than those 
of other structures by more than 20 meV. The distance 
between the BN sheet and the Cu(lll) surface is about 
2.70 A for the Al arrangement and 2.79 A for the A2 
arrangement, which is smaller than the distance between 
adjacent atomic layers in bulk h-BN (3.35 A) but larger 
than that between a BN sheet and the Ni(lll) surface.— 
Therefore, the strength of interaction between the BN 
sheet and Cu(lll) can e expected to be stronger than 



that between the BN layers and weaker than that with 
the Ni(lll) substrate. In the relaxed structure, the B 
and N atom are found not in the same plane, and the 
N atom which is at the atop Cu site is a slighly higher 
(~ 0.02 A) than B atom which is located at the hollow 
site. This rumpling is also much smaller compared to 
that of BN on Ni(lll) substrate where a rumpling of 
0.18 A was observed experimentally. This is also an in- 
dication of relatively weaker interaction between the BN 
layer and Cu(lll) substrate, compare to Ni(lll). 

To find the most stable structure of graphene on 
BN/Cu(lll), self-consistent calculations were performed 
for all possible geometries between the graphene and the 
above two most stable BN arrangements on Cu(lll). 
Relative to the BN lattice, carbon atoms of the graphene 
can be directly over the B and N atoms, or one carbon 
atom in the hollow site and the other above either a B or 
a N atom. Among the six possible heteroepitaxial struc- 
tures, the two atomic arrangements Al (Fig. la) and A2 
(Fig. lb) are energetically degenerate and are predicted 
to be most stable. In both structures, one graphene sub- 
lattice is over the B atom. The other sub-lattice is on 
the hep hollow site of Cu(lll) in Al while it is on the fee 
hollow site of Cu(lll) in A2. The graphene floats over 
the BN sheet at a height of ~ 3.0 A, slightly smaller than 
inter-layer distance of graphite. This confirms that LDA, 
despite the lack of long-range nonlocal correlations, pro- 
duces reasonable interlayer distances in layered van der 
Waals crystal, such as graphite and h-BN, for a deli- 
cate error cancelation between exchange and correlation 
which underlies this apparent performance of LDA.— 
The distance between BN sheet and Cu(lll) changed 
by —0.3 A after adsorption of graphene. However, the 
distance between Cu layers remains essentially the same 
as that in a clean Cu(lll) surface. Further structural 
details are given in Table I. 

Figs. 2(b) and 2(c) show the band structure of opti- 
mized atomic arrangement A2 before and after graphene 
adsorption, respectively. For reference, the band struc- 
ture of clean Cu(lll) surface is shown in Fig. 2(a). The 
rf-bands of copper are apparent in all three band struc- 
tures and are essentially unaffected by BN and graphene 
due to the weak interaction with them, except at the T 
point near the Fermi level. After adsorption of the BN 
sheet, bands characteristic of h-BN appear at about —18 
cV in the band structure while the occupied p-bands of 
BN appear at about 4 eV below the Fermi level at the 
L point. The most interesting feature is the appearance 
of a band in the energy gap of Cu(lll) at K about 2 
eV above the Fermi level. In perfect BN sheet, this va- 
lence band minimum is about 18.4 cV above the low- 
est s-band at K within LDA calculation, but it is re- 
duced to about 18.2 eV after the BN sheet is adsorbed 
on Cu(lll). This shift is due to surface dipole effect of 
Cu(lll) substrate, as indicated by the change of work 
function of the Cu(lll) surface from 5.25 cV to 4.83 eV 
after BN adsorption. Electrons redistribe in response to 
the dipole effect, resulting also to the slight downward 
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shift of some Cu(lll) bands at M and the upward shift 
of some bands at T around the Fermi level. Neverth- 
less, the BN layer remain insulating since the Fermi level 
does not cross any BN bands. After graphene adsorption, 
the Cu(lll) bands at T become completely unoccupied 
and bands characteristic of graphene appear in the en- 
ergy range of —18 cV to —4 eV. The Dirac cone at K is 
located in the energy gap of the substrate, just 0.5 eV 
below the Fermi level. The linear relationship between 
energy and momentum is maintained. This indicates that 
the graphene is n-doped with excess carriers induced by 
the substrate. In addition, a band gap of about 0.2 eV 
opens at K for graphene which is much larger than that 
of graphene on h-BN substrate^, although the interfa- 
cial atomic arrangement between graphene and h-BN is 
the same in both cases. As mentioned above and shown 
in Fig. 1(b), one sublattice of graphene is located directly 
above B or the hep hollow site of Cu(lll) while the other 
is at the hexagonal center of h-BN honeycomb or the fee 
hollow site of Cu(lll). Interaction between the Cu(lll) 
substrate and BN leads to redistribution of charge of BN 
sheet since the BN plane is a little rumpled. Thus the 
gap opening of graphene can be understood based on the 
symmetry breaking the A and B sublattice equivalence, 
which leads to the rehybridization of the valence and 
conduction band states associated with the same Dirac 
point (see Fig. 3(b) and Fig. 3(c)). This gap is essential 
for graphene-based electronic devices for controlling con- 
ductivity. 

Although the Dirac point of graphene is located in the 
energy gap of the Cu(lll) substrate at point K of recipro- 
cal space, the electronic structure of the graphene layer at 
the other points of momentum space is unknown, which 
would have influence on conductivity. The K-resolved 
DOS, which is equivalent to ARUPS experiment, of the 
graphene layer of the A2 structure is shown in Fig. 3(a). 
The band structure is almost the same as that of per- 
fect graphene except the position of the Fermi-level (The 
energy gap cannot be seen clearly at K which is due to 
the resolution in the calculated DOS). No other states 
were observable in this graphene layer. Therefore, the 
electronic property of the graphene layer shows features 
of a typical n-doped graphene. However, it is noted that 
this occurs without an applied gate-voltage^ and the 
n-doped effect is completely induced by the substrate. 
Fig. 3(b) shows the precise K-resolved DOS around K of 
the carbon atom occupying the atop B site and Fig. 3(c) 
shows that of the carbon atom occupying the other sub- 
lattice. Around the K point, the valence band and con- 
duction band of the two sublattices do not touch each 
other due to the symmetry breaking of the two sublat- 
tices. They are different because of the rehybridization 
of the valence band and conduction band states at the 



same Dirac point, as illustrated in Fig. 3(d). For carbon 
over boron, the conduction bands degenerate while for 
carbon in the other sublattice, the valence bands degen- 
erate. The separation between the conduction band peak 
of one carbon and the valence band peak of the other car- 
bon at K is about 0.18 cV which is in good agreement 
with the value determined from calculated band struc- 
ture. The Fermi level is above this gap due to substrate 
induced doping. The valence Dirac cone is populated 
with electrons which have nearly zero effective-mass and 
very high mobility. 

Similar calculations and analysis were carried out for 
structure Al, and its electronic structure is nearly the 
same as that of A2. The only noticablc difference is the 
slightly small band gap, 0.15 cV, which is due to the 
different small difference in atomic structures. 

This n-doped graphene on the BN/Cu(lll) substrate 
can be very useful in realization of graphene-based elec- 
tronic devices. It was demonstrated recently that p- 
type graphene can be easily realized experimentally by 
surface modification with electron acceptor molecule^ 
But a simple and effective method to fabricate p-n junc- 
tion based on graphene is still lacking. Based on re- 
sults of present study, it would be possible to massive 
produce p-n junctions by depositing electron acceptor 
molecules on certain regions of graphene eptaxially grown 
on BN/Cu(lll) substrate. Furthermore, because the gap 
is just 0.5 eV below the Fermi-level, the Fermi-level can 
be easily shifted inside the gap by adsorption of molecule 
or applying a gate voltage to make graphene a real semi- 
conductor. With a 0.2 eV band gap, graphene would 
behave like a narrow-gap semiconductor, such as PbTe 
or PbSe. It may be used as a material in infrared light 
emitting or detecting devices. One example would be 
light emitting quantum dots, which have much more ad- 
vantages than conventional devices. 

In conclusion, we have carried out a first principles in- 
vestigation on graphene grown on Cu(lll) surface with 
a h-BN buffer layer. Electrons transfer from the cop- 
per substrate to graphene through the BN buffer layer 
is predicted which results in a n-doped graphene in the 
absense of a gate voltage. A gap opening at the Dirac 
point just 0.5 eV below the Fermi-level is observed. The 
size of the gap (0.2 eV) is comparable to that of a typ- 
ical narrow gap semiconductor. The Fermi-level can be 
easily shifted inside this gap to make graphene a semi- 
conductor. Considering that the gap size of graphene on 
h-BN or clean Cu(lll) is only a few meV^ it would be 
possible to engineering the gap of graphene through the 
tuning the interaction between the buffer layer and the 
metal surface by choicing different metal surface or buffer 
layer. 
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TABLE I: Interlayer distances, as indicated in Fig. 1, of fully 
relaxed structures. All values are given in A. 
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FIG. 1: The two most stable structures of graphene grown on 
Cu(lll) with a h-BN buffer layer between them. A top view 
and a side view are presented in each case. The Cu atoms 
are represented using big grey spheres and their darkness in- 
creases with the distance from the surface. The carbon atoms 
are represented by small dark spheres. Nitrogen and boron 
atom are represented by medium sized blue and pink spheres 
respectively, do, d\, d 2 and d represent interlayer spacings 
and their values are given in Table I. 



FIG. 2: Band structure of (a) clean Cu(lll); (b) h-BN on 
Cu(lll); and (c) graphene on the BN/Cu(lll). The Fermi 
level is indicated by dash line. 



5 



FIG. 3: (a) K-resolved DOS of the graphene layer, (b) K- 
resolved DOS of carbon atop B along the path indicated in 
(d). (c) Same as (b) but for carbon of the other sublattice. 
(d) Schematic diagram to show the band opening at Dirac 
point due to symmetry breaking between the two sublattice 
(A & B). 



